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^ '. Abstract 
(D ; 

\ We present semi-empirical tight binding calculations on thienylenevinylene oligomers up to the 

c/3 ' hexadecamer stage (n = 16) and ah initio calculations based on the local density approximation 

S ■ up to n = 8. The results correctly describe the experimental variations of the gap versus size, the 

! optical spectra and the electrochemical redox potentials. We propose a simple model to deduce 

I from the band structure of the polymer chain the electronic states of the oligomers close to the gap. 

^ ' We analyze the evolution of the gap as a function of the torsion angle between consecutive cells: 

a^ : 

. the modifications are found to be small up to a ~30° angle. We show that these oligomers possess 

'nI" \ extensive vr-electron delocalization along the molecular backbone which makes them interesting for 

O ■ 

' future electronic applications such as molecular wires. 
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I. INTRODUCTION 



In the microelectronic road map, one approach for the conception of nanodevices is the 
"bottom- up" one which consists of arranging individual atoms or molecules to get a given 
function. This kind of approach has been strongly developed thanks to the scanning tunnel- 
ing microscope which offers the possibility of both observation and manipulation of single 
atoms or molecules (for a recent review, see^). Among the infinite possibilities of synthetic 
organic chemistry, linearly vr-conjugated systems with multinanometer dimensions offer very 
exciting opportunities due to their potential use as molecular wires in nanoelectronic devices. 

Such a purpose requires stable linear 7r-conjugated systems combining optimal vr-electron 
delocalization with approaching the present limit of nanolithographic techniques, i.e. the 100 
A regime. In this context extended vr-conjugated oligomers based on phenyleneethynylene 
and thiopheneethynylene have been proposed as potential molecular wire&-:^. However, all 
these structures exhibit a more or less rapid saturation of the effective conjugation length and 
hence of the HOMO-LUMO gap due to excessive vr-electron confinement related to structural 
factors such as aromatic resonance energy or as rotational disorder-. Consequently, the width 
of the gap converge towards a limiting value for chain lengths significantly shorter than a 
maximum chain dimension. 

Although the first members of oligomers of the thienylenevinylene oligomers (nTVs) have 
been known for several years^'-, their lack of solubility did not allow neither the synthesis of 
longer oligomers nor the detailed analysis of their electronic and electrochemical properties. 
As shown recently, the introduction of solubilizing alkyl chains at the 3 or 3 and 4 positions 
of the thiophene allows the synthesis of highly extended nTVs (up to the hexadecamer stage, 
n = 16) with chain dimensions close to the present limits of the nanolithography (~ 10 nm)^. 
Owing to their recent emergence as a new class of extensively vr-conjugated oligomers, nTVs 
have not been subject to theoretical works, in sharp contrast with parent series of conjugated 
oligomers such as oligothiophenes, oligophenylenevinylenes or oligoanilinies^"— . 

This work presents a theoretical analysis of the chain length dependence of the electronic 
properties of nTVs involving one to 16 thiophene rings. We discuss the nature of the states, 
the optical properties and the ionization energies of the oligomers and the theoretical results 
are then compared to the available experimental data. 
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II. CALCULATIONS 



The atomic and electronic structure of molecules or solids can be calculated by various 
techniques which rely on different degrees of approximation. Ab initio configuration inter- 
action techniques are extremely powerful to predict spectroscopic properties but cannot be 
applied to molecules containing a large number of atoms like nTVs. Thus we have performed 
ab initio calculations based on the local density approximation (LDA)^^ using the DMOL 
code^. For the computation, we use a double numerical basis set (two atomic orbitals for 
each occupied orbital in the free atom) together with polarization functions {2p for H and 3d 
for S and C). All the results presented in this paper are based on the spin-density functional 
of Vosko et ai^^. Results obtained with another functional^ show only minor differences 
(~ 2 meV on the gap of thiophene). To save computational workspace, the hexyl groups 
are replaced by methyl groups ([1]). We check that this does not change the nature and the 
energy of the electronic states of interest, i.e. close to the gap. For thiophene, 2TV and 
3TV, an optimization of the geometry is performed by minimization of the total energy 
with respect to the atomic coordinates-^. We obtain that the molecules are almost planar 
(this point is discussed in the next section) and that the geometry of the different groups 
are almost identical. Thus, for longer nTVs {n > 3), we have just replicated identical cells 
to build planar molecules without optimization of the geometry. Because of computational 
limits, we are able to apply LDA calculations only up to n = 8. For longer chains, a simpler 
approach is required. This is achieved by using a self-consistent tight binding calculation 
which is a semi-empirical technique close to the simple Hiickel theory. C and S atoms are 
represented by one s and three p atomic orbitals, and H atoms by one s orbital. The el- 
ements of the Hamiltonian matrix, restricted to first nearest neighbor interactions and to 
two center integrals, are written in terms of a small number of parameters (15) following 
the rules of ref.— . First nearest neighbor interactions depend on the interatomic distance 
following Harrisons rules^^. For the sake of consistency with LDA calculations, we apply the 
following procedure to determine the parameters. First, we calculate in LDA the electronic 
structure of simple molecules: methane, ethane, ethylene, acethylene and CS2. Then the 
parameters are determined by a fit of the LDA electronic structure using a least square 
minimization technique. The charge transfers between the atoms are also adjusted on those 
obtained in LDA using a MuUiken analysis^^i^. The calculation is made self-consistent by 
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adding to the diagonal matrix elements charge dependent Coulomb terms^. Details on the 
calculations and on the tight binding parameters are given in Appendix. 

Let us discuss the respective limits of LDA and tight binding calculations. LDA tech- 
niques are very powerful to predict ground state properties of molecules or solids^^ but, 
strictly speaking, are not applicable to excited states. Actually, one-electron gaps calcu- 
lated in LDA are usually too small compared to experimental optical gap9i^>22. Thus a 
correction is usually applied to LDA gaps corresponding to a rigid shift of the unoccupied 
states with respect to occupied ones (the so-called scissor operator). The correction results 
from a change of the exchange-correlation potential through the gap (correlation effects) 
and also, for finite systems like molecules, from electrostatic effects and Coulomb interac- 
tion between the electron and the hole in the excited state^^i^. It is often assumed that this 
correction does not depend too much on the size of the chain, i.e. that the variations of the 
gap are well predicted in LDA—. Recent calculations on the gap of silicon nanocrystallites 
tend to confirm this statement^. Finally, the same underestimation of the gap is expected 
here in tight binding because the parameters are adjusted on the electronic structure of 
small molecules calculated in LDA. 



III. RESULTS AND DISCUSSION 



The synthesis and characterization of the optical and electrochemical properties of the 
oligomers have been described elsewhere^. Fig. [2] reproduces typical UV- visible absorption 
spectra recorded in CH2CI2 for 4TV, 6TV, 8TV, 12TV and 16TV. In each case, the spectrum 
exhibits a well-resolved vibronic fine structure indicative of the persistence of a relatively 
planar and rigid geometry even for longest chains. From the spectra we define the optical 
gap by the maximum of the lowest transition which we report on Fig. [3] versus the 
number of sp^ carbons in the molecule {N^ = 6n — 2). Chain extension leads to a shift 
of the maximum absorption corresponding to a narrowing of the gap. It reaches 1.95 eV 
for Nc = 94 {n = 16), the smallest ever reported value for a vr-conjugated oligomer of 
homogeneous chemical structure. A representation of the gap with respect to l/N^ (or 1/n) 
is common in the literature^^ but in our case it is a good approximation only for small 
chains because the gaps saturate for large Nj-. A saturation, corresponding to a limit of the 
effective conjugation length, is reported for all the series of oligomers studied up to now^. 
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Meier et alM have suggested that this hmit is reached when the bathochrom shift (or red 
shift) of Xmax resulting from the lengthening of the chain is smaller than 1 nm. The shift of 
6 nm that we measure between 12TV and 16TV shows that the limit is close to n = 16but 
it is not yet reached. The comparison with other oligomer series shows that nTVs have the 
longest limit of effective conjugation as well as the smallest gap^. It remains to compare the 
variations of the gap with those predicted by the theory. 

Gaps calculated in LDA and in tight binding are plotted in Fig. |3l The results obtained 
by the two techniques are reasonably close which is remarkable considering the simplicity 
of the tight binding calculation. As discussed above, we need to take into account the 
underestimation of the gap by the LDA. Applying to the theoretical gaps an upward shift 
of 1.2 eV (Fig. E]), we obtain a correct agreement between theory and experiments (within 
~0.2 eV): thus the experimental variations of the gap are well described by the theory. 1.2 
eV is a reasonable and common value for the gap correction^*^. The tight binding gap is well 
approximated in eV by 1.85 + j^rn including the shift of 1.2 eV. The limit (1.85 eV) is the 
band gap energy of the polymer (polythienylenevinylene) calculated in tight binding (+ 1.2 
eV), in good agreement with the experimental values of the order of 1.80 eV^^. An exponent 
(1.12) between 1 and 2 is commonly obtained in semiconductor quantum dots^^. 

Fig. m shows the densities of LDA eigenstates from thiophene to 6TV. Similar results are 
obtained in tight binding. We see that the evolution of the HOMO and LUMO states is 
almost symmetrical which implies an efficient delocalization of both states. This is confirmed 
by the analysis of the molecular orbitals showing that these electronic states mainly derive 
from the interactions between vr orbitals along the whole molecule. These results show that 
in principle efficient electron or hole conductivity could be achieved in nTVs providing that 
the n or p-tjpe doping is realized (but to conclude the effect of the charge carriers on the 
conformation of the molecule should be analyzed). 

Fig. [5] shows the electronic band structure of the polymer chain calculated in tight 
binding. The highest occupied band and the lowest unoccupied one have a large dispersion 
in energy which is consistent with an efficient delocalization of the states. Interestingly, 
the energies of these conduction and valence bands are reasonably well represented by 
simple cosine laws: 

= 3.65 - 0.89 cos(A;a/2) (1) 
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eh = 1.13 + 0.98 cos(A;a/2) 



(2) 



where k is the wave vector and a is the length of the translational unit cell which contains 
two thienylenevinylene groups (the energies are given in eV and they do not include the shift 
of 1.2 eV). The formula (1) and (2) are equivalent to the tight binding dispersion laws of a 
linear chain with one orbital in a unit cell of length a/2 (corresponding here to one orbital in 
each thienylenevinylene group) and with a constant coupling only between nearest neighbor 
orbitals. From this we can deduce a simple way to calculate the electronic states of the 
oligomers close to the gap. The main assumption is that the two bands are independent. 
Then the electronic states of nTV are modeled by the eigenstates of the Hamiltonian of 
the finite linear chain containing n orbitals. This is equivalent to the well-known analytical 
Hiickel theory for polyenes^I. The solutions of such problem being analytic, we obtain 
immediately: 



where Cp (hp) represents the electron (hole) states and p is an integer index such that 
1 < p < n . Some electron and hole states {p = 1,2,3) of 6TV calculated with eqn. ([3]) 
and dl]) are represented in Fig. [5] (horizontal dashed lines): they agree relatively well with 
those obtained with the full tight binding calculation (box on the right side of Fig. \5^. The 
agreement is mainly qualitative because many effects are left out in the model such as the 
coupling between the two bands. But this kind of simple model can be particularly helpful 
to investigate different problems related to (like optical properties discussed below). 

We have also calculated the optical absorption of nTVs in tight binding. Working within 
the so-called dipole approximation, the intensity of the absorption is proportional to the 
momentum matrix elements between initial and final electronic states^^. The optical matrix 
elements are computed like in ref.™. Phonon-assisted transitions are not considered in 
our calculations. Fig. |2] shows that the calculated transitions agree relatively well with the 
experiments: in particular we can explain the origin of the two broad bands in the absorption 
spectra (however the energy of the upper band seems to be slightly underestimated by the 




(3) 




(4) 



6 



calculations for large n. Note that this agreement is obtained when a correction of 1.2 eV is 
applied to the gap. Each absorption band is indeed composed of several transitions whose 
number increases with n. The most intense transitions, identified by labels in H] and by 
arrows in [S] in the particular case of 6TV, can be easily interpreted using the results of 
the analytical Hiickel theory for polyenes^^.. The lowest absorption band is composed of 
transitions hp — )■ Cp. The transitions hp — )■ Cm with p ^ m are not visible in Fig. |2) they 
are almost dipole forbidden for the following reasons. The eigenfunctions > and > 
corresponding respectively to the states hp and e^ can be written: 

li^i = E > (5) 

i=l,n 

\rm=Y.<\^i> (6) 

i=l,n 

where i is the index of the thienylenevinylene group. |0f > and |0f > are the effective 
orbitals on the site i corresponding respectively to the valence (h) and conduction (e) bands. 
It is important to point out that the coefficients {a\ or a^) in eqn. (5) and (6) do not depend 
on the band (e or K) because, in the case of a linear chain of identical orbitals with a constant 
coupling between the nearest neighbors, the eigenstates do not depend on the strength of 
the coupling. To characterize the transition /ip — )■ e^, we must write the matrix element of 
the momentum: 

<KIp1C>= E afar < > • (7) 

i j=l,n 

The main contributions in eqn. ([7]) come from the momentum matrix elements between 
orbitals cj)^ > and 0^ localized on the same site, i.e. when i = j. Using the fact that 
A =< \(pi\p\\4>t > is independent of the site i by translational symmetry, we obtain: 

<l^;b1C>= E af ar < 0f . (8) 

i,j=l,n 

where 6pm = 1 if p = m, otherwise. The last equality in eqn. ([8])) which results from the 
orthonormality of the eigenstates explains the selection rule p = m. The same conclusion 
could be obtained using the effective mass theory which has been extensively used in the 
field of semiconductor heterostructures^. The second absorption band (Fig. [1]) is associated 
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with transitions from the states hp to states (e*) composed of s and Pa orbitals locahzed 
on each thiophene group. These latter states have a constant energy which means that 
there is almost no coupling between states localized on neighbor groups. They give rise 
to a flat band in the case of the polymer (Fig. E]). A more detailed comparison between 
theory and experiments is difficult because individual transitions cannot be identified in the 
experimental absorption spectra (Fig. [1]). The broadening of the peaks may result from 
the coupling of the electrons to vibration modes. This interpretation is supported by the 
presence of several equidistant replica of the lowest transition^. The separation of 0.15-0.2 
eV could be associated to stretching modes of the double bonds in the conjugated system. 
Other sources of broadening could be the interactions of the molecule with the solvent as 
well as some rotational disorder. 

In order to investigate the possible effect of the rotational disorder on the gap, we have 
calculated the electronic structure of 2TV when a thiophene group is twisted with respect 
to the remaining of the molecule. The torsion angle is defined in Fig. [H the zero angle 
corresponding to a planar conformation. We have already mentioned that the optimized ge- 
ometry of 2TV corresponds to an almost planar conformation. Because the torsion potential 
is very fiat between -10° and +10° with variations in energy of only few meV a precise deter- 
mination of the torsion angle is clearly beyond the accuracy of such calculation. However, 
in spite of this uncertainty on the position of the energy minimum, our calculations support 
that nTVs are a favorable case to get a planar conformation which allows a maximum delo- 
calization of the electron and hole states. It is interesting to compare with other conjugated 
oligomers. In the case of bithiophene, a ~ 30° angle between the two rings is calculated 
from ab initio calculations^ and ~ 26° from semi-empirical restricted Hartree-Fock Austin 
Model 1 (AMl)^. High-torsion angles are also predicted in the case of three base forms 
of polyaniline, i.e. leocoemeraldine, emeraldine, and pernigranilinei^. A planar conforma- 
tion has been assumed in the case of AMI calculations on oligo(phenylenevinylene)s^, in 
agreement with some fiuorescence studies on stilbene^^ but in disagreement with electron 
diffraction experiments which give an angle of ~30°— . 

The increase in the total energy of 2TV is 0.33 eV at 90° but is only 46 meV at 30°. Thus, 
at room temperature, the thermal fluctuations can probably induce angle variations of the 
order of 20°. Obviously the geometry is also controlled by the interactions of the molecule 
with its environment: for example, planar conformation of oligomers is usually favored in 
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the solid state because of packing effects^!. Fig. |6] shows that the gap energy increases with 
the torsion angle (tight binding calculations predict a smaller increase than LDA). This is 
the expected result since the coupling between vr orbitals of opposite groups is decreasing. 
A similar evolution of the gap has been obtained for biphenyl, bipyrrole and bithiophene^. 
However, the influence of the angle fluctuations on the optical properties is small: a torsion 
angle of 30° gives a shift of only ~ 0.2 eV of the gap of 2TV (Fig. [6]). 

Cyclic Voltammetric experiments have also been performed and analyzed in ref.-. The 
oxidation and reduction potentials of various nTVs (n = 4to 12) are plotted versus 1/Nc on 
(Fig. [7]). We compare these potentials to the ionization potentials and the electron affinities 
calculated for the free molecules. They are obtained in LDA by computing the total energy 
of the molecules in the +1, 0, -1 charge states. The calculated ionization potential of 
thiophene is in good agreement with the experimental one measured in the gas-phase^. 
The experimental oxidation and reduction potentials are given with respect to the Ag/ AgCl 
reference. In order to make quantitative comparisons to free-molecule values, we follow the 
procedure commonly adopted in the literature^^i^: we take the zero of the Ag/AgCl scale to 
be 4.7 eV with respect to vacuum. The calculated ionization potentials and electron affinities 
are approximately linear functions of l/Nf.. This behavior mainly reflects the variation of 
the Coulomb energy associated with the ionization of the molecule. The measured oxidation 
and reduction potentials have a smaller dependence on l/Nc- We attribute the difference 
between experimental and theoretical potentials to solvation effects. To check this, we have 
performed LDA calculations including the Conductor-like Screening Model (COSMO)^'^. 
COSMO is a continuum solvation model where the solute molecule forms a cavity within 
the dielectric continuum of permittivity e that represents the solvent. It provides a good 
approximation of the electrostatic contribution to the solvation energy^^i^. Fig. [7| shows 
that theory and experiments agree well when solvation effects are included. 

In conclusion, we have presented theoretical studies on thienylenevinylene oligomers. 
LDA and tight binding calculations gives results in agreement with the experiments for 
the gap versus size, for the optical properties and for the electrochemical redox potentials. 
We propose a simple model of the electronic structure and of the optical spectra of the 
oligomers which can be particularly useful for further studies. We predict a small effect of 
the rotational disorder on the electronic structure and we obtain an efficient delocalization 
of the electronic states close to the gap. Combined theoretical and experimental results 
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support that nTVs are good candidates to be used as molecular wires. 
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IV. APPENDIX 

The tight binding eigenf unctions tpk of energy are written as linear combinations of 
atomic orbitals i (i is the atomic index and a is the orbital index s, x, y or z): 

^k = Y."'iaXia. (9) 

i,a 

From this we calculate the number of electrons Mj on the ith atom as: 

M = El«Lr, (10) 

k,a 

where the sum over k is restricted to occupied states (including spin). The atomic popu- 
lation Mj are in general different from the neutral-free-atom values Zj, i.e. each atom bears 
a net charge. The self-consistency is incorporated in a simple manner following ref.'iS. We 
assume that only the diagonal terms of the Hamiltonian matrix are charge dependent. They 
are written^^ (in atomic units): 

H.a,a = Hl,^^ + J2{M,-Z,H, (11) 

j 

with 

7., = (4 + U-')-'^' (12) 

Rjj is the distance between atoms i and j. U is the intra-atomic Coulomb energy and the 
quantities H^^ define the s and p atomic levels in the molecule. The nondiagonal matrix 
elements are restricted to first nearest neighbor interactions and to two-centers integrals. 
They are written in terms of a small number of parameters following the rules of Slater and 
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Koster^^. These parameters, fitted on the LDA electronic structure of small molecules, are 
given in[Tl Following ref.— they depend on the interatomic distance d as {d^/dY where do 
is also given in Table 1. The intra-atomic Coulomb energy U is equal to 10 eV for all the 
atoms. 
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TABLE I: Tight-binding parameters following the notations of Slater and Koster ■ . All the energies 
are in electron-volts. The first nearest-neighbor interactions are given for an interatomic distance 
do. 

A-:-/.: C-C C-H C-S 

levels interactions interactions interactions 
Carbon (sso-)=-2.85 (ssa)=-5.50 (sscr)=-3.78 
E,=-14.57 
Ep=-5.67 

Sulfur (spct)=2.93 (spcT=6.90) (spc7)=4.80 
E,=-16.77 
Ep=-9.07 

Hydrogen (ppcr)=3.70 (ppc7)=4.60 
E,=-2.07 

Sulfur (pp7r)=-1.90 (pp7r)=-2.10 
do=1.54 A do=1.07 A do=1.54 A 
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FIG. 1: Molecular structure of 2TV and definition of the dihedral angle investigated in this work. 

FIG. 2: Experiments: UV-visible absorption spectra of various nTVs recorded in CII2CI2. The- 
ory: optical transition energies calculated in tight binding (arrows). The height of the arrows is 
proportional to the oscillator strengths. 

FIG. 3: Experiments: optical gaps of nTVs in solution (+) versus the number of sp'^ carbons 
Nc = 6n — 2. The optical gap of thiophene in solution is from ref.^^. Theory: gaps calculated in 
LDA (o) and in tight binding (x). The calculated gaps are shifted 1.2 eV upward compared to 
experimental data (scale on the right). The line is a fit of the tight binding gaps using 1.85 + jfi^u 
(in eV). 

FIG. 4: Electronic states calculated in LDA from thiophene to 6TV. The gap is represented by 
arrows. 

FIG. 5: Left: band structure of polythienylenevinylene calculated in tight binding (the conduction 
bands are not shifted). Dashed curves: fits of the highest occupied band and the lowest unoccupied 
band by a cosine law. The horizontal dashed lines are the electron and hole states calculated for 
6TV using the simple model of a "s" band. The vertical arrows correspond to the main optical 
transitions. Right: some electron and hole states of 6TV. 

FIG. 6: Gap (continuous line: LDA, dashed line: tight binding) of 2TVS with respect to the torsion 
angle. 

FIG. 7: Electron afhnities (lower +'s) and ionization potentials (upper +'s) of nTVs in the gas 
phase calculated in LDA versus 1/Nc where Nc is the number of sp'^ carbons in the molecule; same 
quantities calculated for the molecules in the solvent (x); experimental reduction potentials (lower 
o's) and oxidation potentials (upper o's). Experimental ionization energy of thiophene in the gas 
phase^(black square). Solid lines: linear fit of the gas phase electron affinities and ionization 
potentials versus l/Nc- 
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